Mott glass in site-diluted S — 1 antiferromagnets with single-ion anisotropy 
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The interplay between site dilution and quantum fluctuations in S=l Heisenberg antiferromagnets 
on the square lattice is investigated using quantum Monte Carlo simulations. Quantum fluctuations 
are tuned by a single-ion anisotropy D. In the clean limit, a sufficiently large D > D c — 5.65(2) J 
forces each spin into its ms = state, and thus destabilizes antiferromagnetic order. In the presence 
of site dilution, quantum fluctuations are found to destroy Neel order before the percolation threshold 
of the lattice is reached, if D exceeds a critical value D* = 2.3(2) J. This mechanism opens up an 
extended quantum-disordered Mott glass phase on the percolated lattice, characterized by a gapless 
spectrum and vanishing uniform susceptibility. 

PACS numbers: 75.10.Jm, 75.10.Nr, 75.40.Cx, 64.60.Ak 
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Strongly interacting quantum systems on random lat- 
tices offer the possibility of realizing genuine quantum 
phases with unconventional properties lj. The inhomo- 
geneous nature of such disordered lattices can lead to the 
existence of strongly correlated regions with excitations 
of arbitrarily low energy, and yet an overall absence of 
long range order. Therefore quantum "glassy" phases 
are possible which, at the same time, have a quantum- 
disordered ground state and a finite low-energy density of 
states, not associated with any Goldstone mode [l|, 12, Ijj . 

A large variety of physical systems potentially lends it- 
self to realize the theoretical scenario of quantum glassy 
phases and quantum phase transitions from such phases 
into more conventional long-range ordered ones. Rele- 
vant examples are 4 He films on disordered substrates £|, 
type-II superconductors with columnar defects []| and 
trapped cold atoms in disordered optical potentials 
Unfortunately, only in few cases has possible evidence for 
quantum glassy phases been observed, and the agreement 
with theoretical predictions is still unsatisfactory. 

In this paper we numerically demonstrate the emer- 
gence of a highly non-trivial interplay between lat- 
tice randomness and quantum fluctuations in low- 
dimensional antiferromagnets. Disorder can typically be 
introduced in antiferromagnets by dilution of the mag- 
netic lattice with non-magnetic ions, and an extremely 
high level of control on the concentration of dopants can 
be reached Q. Recent experiments on site diluted com- 
pounds with antiferromagnetic interactions have shown 
that the destruction of magnetic long-range order (LRO) 
due to doping occurs at the percolation threshold of 
the lattice 0, I n t ne case of the S — 1/2 square- 
lattice Heisenberg antiferromagnet this is clearly under- 
stood since this particular system is too far from a quan- 
tum critical point to develop quantum critical fluctua- 
tions triggered by site dilution 10] . In this paper, we 
consider a realistic magnetic model on the square lat- 
tice in which the strength of quantum fluctuations can 
be tuned at will, driving the system from a classically 
ordered state to a quantum disordered one. Non-linear 



quantum fluctuations enhanced by lattice disorder give 
rise to a novel quantum disordered phase which is gap- 
less and has a vanishing uniform susceptibility. An exact 
bosonic mapping of the model allows us to identify this 
regime with a Mott-glass phase [j| of a commensurately 
filled lattice gas of interacting bosons. The application 
of a magnetic field to such a phase drives it into a disor- 
dered Bose-glass phase 0, which maintains the gapless 
spectrum but acquires a finite susceptibility. 

We investigate the S=l square-lattice anisotropic 
Heisenberg antiferromagnet (SLAHAF) with site dilu- 
tion, whose Hamiltonian reads 

H .iy^< l( ,S,-S, +£^e l (Sf) 2 -/ l J^e l Sf. (1) 

(ij) 1 i 

Here S denotes S—l spin operators, and (ij) enumer- 
ates pairs of nearest neighboring sites on the square lat- 
tice. The variables ej take the values or 1 with probabil- 
ity p and 1 — p respectively, where p is the concentration 
of non-magnetic dopants. We make use of the Stochas- 
tic Series Expansion Quantum Monte Carlo (SSE-QMC) 
method based on the operator-loop algorithm jjjj , which 
allows us to faithfully monitor the T — physics by a (3- 
doubling approach yj on L x L lattices with L up to 36 
sites. The results are typically averaged over at least 200 
disorder realizations. 

In the clean limit p = 0, and at zero magnetic field 
h = 0, the SU(2) symmetric version of this model for 
D = is known exactly to show Neel LRO 12] . The 
presence of finite single-ion anisotropy (D > 0) reduces 
the symmetry to U(l), and antiferromagnetic ordering 
takes place in the xy plane. In the limit D/J 1 the 
ground state becomes |^o) = Yli \ m S — 0}i to minimize 
the anisotropy term. Such a state has no antiferromag- 
netic order, since all spin-spin correlation functions are 
simply zero. Therefore, at a critical ratio (D / J) c a quan- 
tum phase transition |13| occurs between the xy-ordered 
regime and a gapped spin-liquid state with short-range 
correlations. Our QMC simulations of Eq. with h,p — 
provide an estimate of (D/J) c — 5.65(2) J, obtained by 
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linear scaling of the correlation length for the x(y) spin 
components ^ a:x (w) ^ £ The quantum-critical scaling 
of the static structure factor S xx ^ yy '(TT,Tr) ~ L 1 I V ^ Z at 
the above point satisfies the expected 3d XY universality 
class with 7 « 1.32, v w 0.67 and z = 1 Q. 

A better insight into the nature of the quantum phase 
transition tuned by the anisotropy is achieved by exactly 
mapping the spins onto bosons through the Holstein- 
Primakoff (HP) transformation S± = y/l — th/2 bi and 
Sf = 1 — rii (together with a ir rotation of one of the two 
sublattices), which gives 



- ]T J« (« t - l)( nj - 1) + D J] ei(ni - l) 2 
-hJ^^CiHi. (2) 
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FIG. 1: (color online) Phase diagram of the site-diluted 
5=1 SLHAF. The magnetic phases (and the corresponding 
bosonic ones) are indicated. QD = quantum disordered. The 
mean-field (MF) line corresponds to the low-dilution limit 
(D/JJcj, = (D/J) c , (1-p). 



Here rij = bjbi — 0, 1, 2 is the dynamically constrained 



occupation number, and Jj. 



JciCj prevents bosons 



from hopping onto doped sites and from experiencing re- 
pulsion from those sites. It is evident that, apart from 
the square-root terms in the hopping Hamiltonian, the 
model of Eq. J5J is a Bose-Hubbard model with soft-core 
interactions, allowing for up to 2 particles per site. For 
h = the Z2 symmetry of the spin model translates into 
a particle-hole symmetry of the bosonic model, which im- 
plies that the system is exactly half- filled, (rij) = 1. Inter- 
estingly, this still holds true in the presence of site dilu- 
tion, which does not alter the particle-hole symmetry of 
the Hamiltonian. In particular, the anisotropy term D of 
the spin model translates into an on-site repulsion for the 
HP bosons, such that the transition driven by increasing 
D/J can be understood as a superfluid-to-Mott-insulator 
(SF-MI) transition driven by the ratio between repulsion 
and hopping 0. The observed (d + 1)-XY universality 
class is consistent with the commensurate filling of the 
lattice [l|. 

Introducing site dilution into the Hamiltonian Eq. l(T)l 
hence offers the intriguing opportunity of studying disor- 
der effects on a SF-MI transition, both at commensurate 
(ft = 0) and at incommensurate fillings (ft ^ 0). In this 
paper, we mainly focus on the ft — regime, and we 
will discuss the case of ft 7^ in a forthcoming publica- 
tion Moreoever, the model of Eq. QJ lends itself to 
studying the behavior of antiferromagnetic order around 
the percolation threshold of the site-diluted square lat- 
tice p = p* = 0.407253... |l9j under continuous tuning 
of the strength of quantum fluctuations. In the limit 
D = 0, it has been demonstrated numerically that the 
5=1 SLHAF on a site-diluted lattice retains LRO up to 
the percolation threshold . The system at finite p and 
finite D interpolates then between the 3<i-XY transition 
at p = and the percolation transition at D = 0. 



In Fig. ^ we show the complete phase diagram of 
the system in the p — D plane at T = and ft = 0. The 
boundary lines are estimated as above, using the criterion 
£xx{yy) ^ ^ This phase diagram suggests that the in- 
terplay between disorder and quantum fluctuations leads 
to a significant departure both from the SF-MI transi- 
tion at p = and from the percolation driven transi- 
tion at D = 0. Let us start from the limit D = 0, and 
focus on the behavior of the system at the percolation 
threshold p = p* . Here we clearly observe that, upon in- 
creasing the ratio D/J, the transition line from magnetic 
LRO to disorder passes through a multi-critical point at 
(D/J)* = 2.3(2), beyond which it departs from the per- 
colation threshold p = p* and bends towards systemati- 
cally smaller p values for larger D/J. Starting from the 
opposite limit p = 0, we observe that an infinitesimal 
amount of disorder leads to a shift to lower values of the 
critical anisotropy (D/ J) c that causes the destruction of 
LRO. Putting together these two pieces of information, 
we can conclude that for any finite disorder concentra- 
tion p there is a non-trivial disorder-dependent value of 
the anisotropy (D/J) CjP < (Z?/J) c ,o at which non-linear 
quantum fluctuations, enhanced by the reduced connec- 
tivity of the diluted lattice, are able to destroy the LRO. 

The interplay between disorder and quantum fluctu- 
ations can be quantitatively understood in the relevant 
limit of p — p* . At this point, the percolating nature 
of the largest cluster in the system crucially depends 
on Id links connecti ng q uasi-2(i islands with higher local 
connectivity (blobs) |2l|. For D/J, well below (D/J) c ^ 
the quasi-2<i blobs have a well-defined staggered moment 
(see Fig. which fluctuates with a characteristic fre- 
quency given by the finite-size gap of the blob if the 
blobs are disconnected. When the blobs are connected 
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FIG. 2: (color online) Ground-state correlation length of the 
5=1 anisotropic Heisenberg antiferromagnet on the chain as 
a function of the anisotropy. The illustrations sketch the or- 
dered/disordered phase on the percolating cluster when the 
ID correlation length is larger /smaller than the average link 
length (Zo) (see text). The arrows on the 2d blobs indicate the 
local staggered moment ~52 iehloh (—iy(S^ y ^} or, alternatively, 
the local phase of the bosons. 



with their average linear size, analogous to what hap- 
pens at p = p*. For weak dilution p <C p*, on the 
other hand, the disordered lattice can be approximately 
represented as a homogeneous lattice with reduced ef- 
fective coordination z e g = z(l — p), corresponding to 
the Hamiltonian of Eq. Q in which the random vari- 
ables are substituted by their average (e») = 1 — p. In 
this case the critical value (£>/</) CiP should be linearly 
shifted by the disorder with respect to the p = case, 
(D/J) c , p w (1— p) (D/J) c> q, which is in very good agree- 
ment with our data for p < 0.05 (see Fig. 



(a) 10 



a L=28 
• L=36 
v L=40 





5 10 15 20 25 
h 1 



by links to form the percolating cluster, their local stag- 
gered moments are either long-range correlated or not, 
depending; on the decay of correlations along the Id links 
[22j. If the characteristic length Zo of the Id links is 
shorter than the correlation length for the S = 1 
anisotropic Heisenberg model of Eq. Q on a chain, the 
links are able to establish 2d LRO in the system. Oth- 
erwise the blobs are uncorrelated, and the system enters 
a quantum-disordered phase. An earlier detailed study 
\2l\ of the geometry of the percolating cluster reports 
an average link length (Zo) ~ 2.7 for site percolation 
on the square lattice. We have studied the correlation 
properties of the anisotropic Hamiltonian Eq. Q on a 
chain with SSE-QMC, systematically estimating the one- 
dimensional equal-time correlation length as a func- 
tion of the anisotropy D/J. According to the above argu- 
ment, the loss of LRO on the percolating cluster should 
occur for (D/J)* such that £\d[{D/J)*] w (Z ). This cri- 
terion leads to the estimate (D/J)* rj 2.2 (Fig.[5J, which 
is in excellent agreement with the location of the mul- 
ticritical point estimated by QMC, as shown in Fig. 1. 
This means that we quantitatively understand the devia- 
tion of the magnetic transition from percolation in terms 
of a critical enhancement of local quantum fluctuations 
on the weak links of the percolating cluster. 

When p <p* , the links connecting quasi-2<i islands ac- 
quire a higher connectivity, evolving from single chains to 
decorated chains or ladders. The critical value (D/J) c<p 
for p < p* corresponds to the one at which the corre- 
lation length on the quasi-ld links becomes comparable 



FIG. 3: (color online) Magnetic response in the Mott glass 
phase (p = 0.2, D/J — 4.4): (a) uniform susceptibility, (b) 
uniform magnetization upon application of a weak field. In 
the inset: different magnetization curves starting from the 
XY ordered phase (D/J = 3), Mott glass phase (D/J = 4.4), 
and QD gapped phase (D/J = 5.8). For D/J — 3 we have 
plotted m/2, and, for D/J = 5.8, h is rescaled by a factor of 
3. Here L = 28 and p = 0.2. 

The deviation of the magnetic transition from the 
percolation transition opens a novel quantum-disordered 
phase that exists on percolated lattices with connectiv- 
ities arbitrarily close to the clean system as D/J — > 
(D/J) c fi. This phase has remarkable unconventional 
properties, as evidenced by its response to a uniform 
magnetic field. In Fig. OJa) we show the low-temperature 
uniform susceptibility Xu(T —> 0) for a representative 
point p = 0.2, D/J = 4.4. It is observed that Xu(T -> 0) 
tends to 0, as in a phase with a singlet ground state, so 
that, in the bosonic language, this phase is incompress- 
ible at T = 0. Yet the specific temperature dependence 
of Xu appears to follow an unconventional exponential 
law Xu(T) ~ exp[—y/ A/T\ for more than two decades of 
the QMC data. This is clearly different from the conven- 
tional behavior (~ exp[— A/T]) expected in the presence 
of a finite singlet-to-triplet gap A. To probe this gap 
directly, we study the response of the system to small 
but finite uniform fields. Fig. Efb) shows the low-field 
magnetization curve of the system in the quantum dis- 
ordered phase. This is strong numerical evidence that 
the magnetization is finite for arbitrarily small fields, 
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although it grows slowly with h, following the uncon- 
ventional exponential behavior m(h) ~ exp[— A/h] over 
about two decades of the QMC data. The magnetiza- 
tion clearly shows that the singlet-to-triplet gap vanishes 
even at h = 0, and this observation is still compatible 
with a vanishing uniform susceptibility for h = be- 
cause Xu — dm/dh ~ h~ 2 exp[— A/h] — > for h — ► 0, 
but it is non-zero for any finite /i. The magnetization 
curve in the novel disordered phase contrasts with the 
linear behavior of m(h) when starting from the XY or- 
dered phase and with the gapped behavior when starting 
from the Mott insulating phase (inset of Fig. |3J). The 
exponential T-dependence of \u and ^.-dependence of m 
can be fully captured by a simple model in which the re- 
sponse to a weak field is given by independent rare clean 
regions of the percolating cluster, whose exponentially 
rare statistics reflects itself in the exponentially small re- 
sponse functions 0,H3- The gapless nature of the spec- 
trum is a property of the entire novel disordered phase 
for (D/J) c>p < D/J < (D/J) C:0 . In fact, for any p one 
can always find an arbitrarily large (albeit rare) clean 
region which locally approximates the behavior of the 
clean square lattice and which has consequently arbitrar- 
ily small local excitations as long as D/J < (D/J) Ci q. 

To summarize, we have seen evidence of an incompress- 
ible, yet gapless phase of the system. The amorphous 
nature of this disordered phase and its rich low-energy 
dynamics make it akin to glassy phases, and at the same 
time its incompressibility, along with the commensurate 
filling of the lattice, bears a strong resemblance to a Mott 
insulating phase. Hence it can be called a Mott glass, a 
name that has been used in the recent literature for dirty- 
boson models at commensurate filling Q in analogy with 
a phase of disordered fermions having the same thermo- 
dynamic/spectral signatures [2j. Moreover, for any finite 
h the ground-state phase acquires a finite susceptibility, 
which means that the Mott-glass phase is driven into a 
Bose- glass [lj which is disordered, gapless and compress- 
ible Eg. 

The occurrence of a bosonic Mott glass phase in the 
in site-diluted S = 1 antiferromagnet with single-ion 
anisotropy is well understood within the bosonic map- 
ping of Eq. @ 0- In particular we wish to stress 
that the Hamiltonian of Eq. QJ, both in its clean and 
disordered form, is a reliable description of the mag- 
netic degrees of freedom in strongly anisotropic S—l 
antiferromagnetic insulators with and without doping. 
In fact, the gapped quantum-disordered state induced 
by the large anisotropy term is realized in a variety of 
Ni-based compounds, such as Ni(C 5 H 5 NO) 6 (N03)2 O, 
[Ni(C 5 H 5 NO) 6 ](C10 4 )2 HI, and the more recently in- 
vestigated NiCl 2 -4SC(NH 2 ) 2 17]. These systems have in 
general a three-dimensional magnetic lattice, but the re- 
sults we have shown for the 2c? case can be generalized 
straightforwardly to 3d. Moreover one can imagine tun- 



ing the D/J ratio experimentally by applying hydrostatic 
pressure to the crystals. 

The Zi symmetry of the magnetic Hamiltonian in zero 
field for all these systems translates into particle-hole 
symmetry and commensurate (half) filling of the soft- 
core bosons on any lattice site, regardless the geometry 
of the lattice itself, namely also on random percolating 
clusters. This ingredient, which is crucial for the appear- 
ance of the Mott glass, appears hard to be achieved in 
systems of real bosons, as e.g. in bosons in optical lat- 
tices, where the introduction of disorder will invariably 
change the particle occupation near lattice defects. 
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